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1. Introduction 

Materials exhibiting commensurately and incommensurately structurally modulated 
phases were extensively investigated during the last years. BCCD, (CH 3 ) 3 NCH 2 COO • 
CaCl2 • H 2 0, is an outstanding example: it shows a rich variety of phases, and exhibits 
many intricate phenomena in its phase diagram, such as indications for accumulation 
points of structure branchings, and the so-called dielectric Tg-anomaly []J . The problem 
of the modulation and its temperature-dependence still poses many open questions. 
The plenitude of experimental results led to a continuing theoretical interest, and many 
different models for their explanation have been developed (for reviews see references 

ID- 

We recently used the concept of symmetry-adapted local modes || to derive the 
Double Ising Spin (DIS) model [§, |6j as a symmetry-based pseudo spin model with 
two two-component pseudo spin variables per crystallographic unit cell: the symmetry- 
breaking atomic displacements occuring at the transition from the high-temperature 
(high-symmetry) para phase to the modulated phases at lower temperatures are 
expanded in terms of an appropriately chosen basis set, and the pseudo spin components 
r and a are the signs of the relevant coordinates 0, ||. In the case of materials 
with high-temperature phase space group Pnma (e.g. BCCD), this approach yields 
the Hamiltonian of the DIS model 

H = TijkTij^k+l) + aijk(Tij(k+l) + — X! (&ijk T ij(k+l) ~ r ijk<^ij(k+l)j 

ijk ijk ijk 

+ Tijk [j{i+l)jk + T i(j+l)k) + J' X! a Vk {^(i+^jk + 0"i(j+l)fc) • (1) 

ijk ijk 

i, j, k label the pseudo spin positions along the a, b, c-directions respectively. At 
T = 0, the in-layer couplings J < 0, J' < lead to a ferro-ordering of the layers 
perpendicular to c. In c-direction, frustrations and therefore modulations arise because 
of the antagonistic effects of the symmetric nearest neighbour interactions K and L on 
the one hand and the antisymmetric interaction M > on the other hand. The stable 
structures (pseudo spin profiles) obtained from the statistical mechanics treatment of H 
are characterized by their wave numbers q given in units of 2n/ (pseudo spin spacing). 

At T = 0, five stable phases are separated by multiphase lines in the K/M-L/M- 
plane. Along these lines, an infinity of different phases are degenerate. We restrict 
ourselves to ferroelectric r-r-interactions, i.e. K < 0; the phase diagram for K > 
is obtained by exploiting the properties of H . In the following, we will consider the 
multiphase lines (a) and (b) separating, in the notation of reference 0, phase V (the 
pseudo spin structure repeats itself after four layers), from the ferroelectric phase I (all 
r- and cr-spins in ferro-order) and the mixed phase III (r-spins in ferro-, cx-spins in 
antiferro-order along c) respectively. 
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For low temperatures, all phases degenerate at line (a) (with K < and L < 0) 
and obeying some rules have a finite stability range in the vicinity of this line [0]. 
Incidentally, this also holds for the four-state chiral clock model which can be derived 
from the DIS model for the case K = L < and J = J' < || . 

In contrast to this, in the vicinity of line (b) (with K < and L > 0) only four 
different phases (characterized by the wave numbers q = 0, |, |, and |) are stable at 
very low temperatures || . Other phases are created at higher temperatures by structure 
combination branching processes. 

In the present paper, we report on a remarkable type of transitions in the DIS 
model for the case K < and L > 0: these 'internal' transitions occuring at 
certain temperatures are characterized by a structural change in the pseudo spin 
arrangement which is not accompanied by a change of q and - in some cases - not even 
by a change of the symmetries of the profile. 

The paper is organized as follows. In section ^ the pseudo spin profiles calculated in 
mean field approximation (MFA) are discussed. The results are substantiated in section 
[3] by employing the mean field transfer matrix method and Monte Carlo simulations, 
which prove that these transitions are no artefacts of MFA. In section f|, the internal 
transition in the q = |-phase of the DIS model is used to explain a dielectric anomaly in 
the fourfold phase of BCCD as well as seemingly contradictory experimental structure 
analysis results. Section [5] contains a summarizing discussion. 



2. Mean field investigation of the internal transitions 

The mean field approximation is often used for analysing the finite temperature 
behaviour of pseudo spin models. Starting point is the variational Hamiltonian 

H = -J2lL k - S ijk (2) 

ijk 

for systems with a modulation in (OOl)-direction, with = (r^fc, cr^). The product 
V k 'S_ijk describes the interaction of the spin 5^ (which in our case has two components) 
with a site-dependent mean field rj k depending on the averages of the neighbouring spins. 
The mean values ti = (r^) and si = (<Jije) of the pseudo spins r and a in the layer I 
are determined by the MFA equilibrium equations (/3 = j^f): 

tt = tanh ((3r)i t t) , 

si = tanh (Pr]2 t e) , (3) 
with the mean fields 

M 

r} U = - 4Jt £ - K (t t+1 + ti-x) + — (s i+ i - st-i) , 
V24 = ~ 4J's £ - L (s t+1 + St-^) - — (t e +i - ^-l) • 
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Figure [I] shows a typical mean field phase diagram in the plane A = 0.1, j — f = —0.25, 

K \ - A 

Mi Mi 



where we have introduced the reduced interaction parameters k = tj, A = tt, j 



m~i ? ~ If anc ^ K - = — A); 6* = ^fj- denotes the reduced temperature. The shaded 
areas contain incommensurate and higher order commensurate phases. The different 
commensurate phases are characterized by their wave numbers. Para phase, ferro phase 
and modulated phases meet at the Lifshitz point (LP). At 9 = the considered cut 
through the temperature-interaction phase diagram intersects the multiphase line at 
the multiphase point MP. 

The existence of accumulation points of structure branchings in the case K < 
and L > was verified || by formulating the mean field equilibrium equations (3) as a 
four-dimensional mapping and using a fixed-point technique analogous to that employed 
for the ANNNI model [10|. Also calculated were the self- and the interaction-energies 
of discommensurations. 

As a consequence of the signs of the parameters K, L and M (i.e. K < 0, L > 0, 
and M > 0), all phases degenerate at the multiphase line (b) consist of bands of two 
types of two-layer sequences viz. (+ ^) and (^ ~), where the upper (lower) symbols 
represent the signs of the r- (a-) spin. We introduce a bracketed notation {Ai . . . A k ) 
for the different phases: for example, (23) denotes a repetitive pseudo spin structure 
consisting of a band with two (+ ^) sequences followed by a band with three (^ ~) 
sequences. The wave number for the above phase is 



k 

k 

4 E A 

i=i 



The equilibrium pseudo spin profiles for configurations with a fixed wave number q at a 
given temperature are obtained by finding the self-consistent solutions of the equations 
(3). This is done by starting from different initial spin configurations with the same 
modulation period and iterating. The respective free energies are calculated and the 
solution with the lowest free energy describes the stable spin configuration. Keeping 
A, j, and j' fixed, the reduced parameters 6 and k_ are then varied in such a way 
that the resulting line 9 (k_) always lies wholly in the stability region of the phase 
under investigation. Figure shows the obtained spin profiles for the phases (12) (wave 
number q = \) and (2) {q = |) with A = 0.05 and j = j' = — 1. For both wave 
numbers there is a discontinuous internal phase transition. At higher temperatures, the 
q = |-phase exhibits also a continuous transition to a modification with two disordered 
layers. 

The obtained spin profiles for the (2)-phase (for layers 1 to 8, i.e. for one modulation 
period) exhibit symmetries given by 
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(a) at low temperatures: 

t : t\ £2 *2 ti —t% —t2 — £2 — t\ 

s : —S\ S2 — S2 S\ S\ —S2 S2 —S\ 

(b) at high temperatures: 

t : ti £2 ^2 ti ~ti — ^2 —t\ 

s : — Si — S 2 S 2 Si Si S2 —S2 —Si. 

At a reduced temperature 9- m i the cr-spins averages change signs in layers 2, 3, 6 and 
7, whereas the symmetry elements of the spin profiles are left unchanged. Furthermore, 
also the r-spins (in all layers) show discontinuities at 9^. 

For the (12)-phase, the spin profiles in the vicinity of the discontinuous transition 

are: 

(a) at low temperatures: 

t : t\ t\ —t2 —tz —£3 —t2 
s : -St Si s 2 -s 3 s 3 -s 2 

(b) at high temperatures: 

t : ti t 2 £3 —ti —t 2 —£3 
s : -Si s 2 s 3 Si -s 2 -s 3 . 
In this case the r-spins in layer 3 (and the cr-spins in layers 4 and 5) change signs, 
leading to a change of the symmetries of the profile. 

In general not all combinations of the signs of the spins were taken as possible 
starting points for solving the equilibrium equations self-consistently. However, we 
substantiated our results for selected values of the parameters by considering all 2 11 
(phase (12)) or 2 15 (phase (2)) possible sign combinations (changing the signs of all 
spins leaves the free energy unchanged). 

In order to characterize the different modifications, a Fourier analysis of the spatial 
modulation of te and Si was performed. In the following, the results for the q = |-phase 
are discussed, as they will be needed for the application to BCCD in section |j. For 
reduced temperatures slightly below {6 = 2.967) and above {6 = 2.993) the reduced 
temperature 6^ nt of the internal transition, the pseudo spin averages (profiles) for one 
period of the spatial modulation are shown in figure together with their harmonics. 
Due to the symmetries in the pseudo spin profiles, the Fourier components of ti (s^) are 
all odd (even) about i — |, and only odd orders in sine (cosine) occur in the Fourier 
series: 



7T I 3"7r 1 

t e = aisin(-[*- -]) + a 3 sin(— [I- -]) 
s t = h cos(^ - + 63 cos(^[£ - -]). 



(4) 
(5) 
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Upon cooling (see figure [|), the relative strength of the third order harmonic with 
respect to the first order harmonic increases smoothly from zero near the transition 
to the para phase to about one third near 9- m ^ (both for the t- and s-profiles). When 
dropping below 0j nt , only a slight increase in a\ and a 3 but remarkable changes in 61 
and 63 occur: 63 jumps from 0.391 at B = 2.993 to -1.089 at 6 = 2.967 while |&i| 
decreases from 1.136 to 0.536. Below 9[ n i, the s-profile is governed by the strong third 
order harmonic (I&3/&1I > 2.03), which continuously increases to the ground state value 
I&3I = 1.307, whereas the first harmonic keeps its nearly constant value \bi \ = 0.54. 

In MFA, these discontinuous internal phase transitions were obtained for all 
investigated phases of the DIS model, execpting only the (1)- and the (00}- (ferro) 
phases. They occur for all values of A and k as long as A > and k < 0, as well as for 
all in-layer couplings j and f (even for very large values). 



3. Verification by mean field transfer matrix and Monte Carlo methods 

In the previous section, the DIS model was treated in the mean field approximation. 
The question arises whether the observed internal transitions are an artefact of this 
approximation or if they are also present if one goes beyond mean field theory. 

This problem should indeed be considered: in MFA, the Axial Next Nearest 



Neighbour Ising (ANNNI) model [11| with weak in-layer couplings exhibits, for example, 



phase transitions from the low temperature phases to high temperature asymmetric and 



partially disordered phases |L2]. Nakanishi [13| could show that these phases are MFA 
artefacts by demonstrating that they do not occur if the mean field transfer matrix 
(MFTM) method is used, in which the interactions in direction of the modulations are 



treated exactly. The same conclusion was drawn from Monte Carlo simulations Jl4 
In the MFTM method, the variational Hamiltonian for the DIS model is 

Ho = T ijkTij(k+l) + L <Jijk&ij(k+l) 

ijk ijk 

+ Y ^2( a ijknj(k+X) - crij(k+l)Tijk) - S-ijk- (6) 

ijk ijk 

In the MFA variational Hamiltonian (1) all interactions are described by mean fields. 
In contrast, in equation (6) site-dependent mean fields are only introduced for the 
interactions perpendicular to the direction of modulation. As the two-component mean 
fields r/ k do not depend on the indices i and j the problem is reduced to one dimension 
and the indices % and j will be omitted. 

Using the Bogoliubov variational principle, the MFTM free energy is given by 

F = 2,7^4 + 2J' 5>| + £ Vi,e te + Yl W t - k B T In Z (7) 
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with 



Z = Tr exp (k ^ r m + Lj2°e + y ^ ( r m ^ - ^+1 r £) + ( h U T z + k*4 a e) ) ( 8 ) 



and K = L = /3L, M = pM and & M = pr] 1A , h 2/ 
the partition function Z we introduce the transfer matrix 

/ e h^+h^ o \ 

e hi,i~h2, e o 

e -h 1>e +h 2 , e o 





In order to calculate 



Tt 



V o 







-hi. 



( e K+L 
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K-L-M 



-K+L+M 



,K-L+M 
3 K+L 



,-K-L 



-K+L-M 



-K-L 



-K-L 



3 K+L 



V 



K-L 



-K+L-M „K-L+M 



3 -K+L+M 
,K-£-M 



3 K+L 



For a commensurate phase with period p the mean fields repeat themselves after p layers. 

The partition function (R) is then given by Z = where a max is the maximum 

v 

eigenvalue of the matrix product n Tt and N is the number of layers in the crystal. 

i=i 

The mean values t? and se are calculated by 
d In aw-* 1 



t, 



Si 



dhj 
dlna m 



a 



■max 
1 



a,. 



a. 



a, 



a, 



T\...Ti-xH Tg...T p 
Ti...T^ifi s T£...Tp\ 



(9) 



(a max | and 
eigenvalue a r 



&max) are the left and the right eigenvectors corresponding to the 
x , whereas the matrices // and fi s are given by 



The mean fields 
Vi r 



( 1 



V o 





1 











and 



\ 



-1 J 



V2, 



and 



f 1 








\ 





-1 














1 





V o 








-1 J 



(10) 



result from the minimization of the free energy (0). Equations (9) and (|10"D are to be 
solved self-consistently. 

As an example, we will discuss the application of this method to the q — | phase. 
Calculations in mean field approximation (see previous section) yield a first order phase 
transition at 8^ between a low- and a high-temperature modification (figure |2[). These 
two spin profiles are used as starting points for the self-consistent solution of the MFTM 
equations (9) and (|T0|). 

Figure |5| shows the resulting free energies obtained by varying the model parameters 
in the same manner as in the previous section, i.e. one proceeds along the same 
curve as for the MFA calculations. It is clearly seen that the free energy curves of the 
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two modifications intersect at 0^ = 2.91, yielding a first-order phase transition. As 
expected, the temperature 9^ is slightly shifted to lower temperatures as compared to 
the MFA result. 

In the MFTM method some pseudo spin interactions, though not those in 
modulation direction (which should be the most important ones), are still replaced 
by an interaction of the pseudo spins with a mean field. To make sure that this 
approximation is not responsible for the occurrence of internal phase transitions, Monte 
Carlo simulations were performed for phases with a short modulation length. No 
attempt was made to determine the global phase diagram. 

Systems of sizes L x L x M were simulated, with L ranging from 10 to 40. The 
number of layers, M, is given by the modulation length of the phase under investigation, 
e.g. M = 6 for the q = |-phase. This is done in order to avoid the development of 
fluctuations resulting from (meta) stable phases with a larger modulation length. The 
standard one-spin-flip Metropolis algorithm is used. 

Figure || shows the calculated specific heat as a function of the reduced temperature 
for the q = |-phase, for interactions given by A = 0.1, k — —0.57, and j = f = —0.5 and 
a size of 30 x 30 x 6. There are two different peaks. The peak at higher temperatures 
corresponds to a second order phase transition to the disordered phase. In the present 
context, the narrow peak is the interesting one. It corresponds to the internal first order 
phase transition. The resulting spin profiles for temperatures below and above the peak 
position compare favorably with the profiles obtained in MFA. In the inset showing the 
calculated energy the discontinuous character of the internal phase transition is clearly 
visible. 

4. Application to BCCD 

Due to the pseudo symmetry of half a lattice constant along c, BCCD can be 
considered to be built up of successive layers of half cells perpendicular to c, which 
we label by the subscript I. Applying the method described in reference |]] to BCCD, 
the symmetry-breaking atomic displacements occuring below the transition from the 
unmodulated high-temperature phase to the modulated phases are expanded in terms 
of a localized symmetry-adapted basis set; the respective coordinates (mode amplitudes) 
are projected onto two- valued pseudo spin variables leading to symmetry-based pseudo 
spin models. Every basis vector (symmetry-adapted local mode, SALM) describes a 
collective displacement of all atoms in a half cell. One obtains two relevant SALMs 
for the front and two for the back half cell. Adequate superpositions of these describe 
structural modulations transforming according to the irreducible representations A2 and 
A3 of the group of the wave vector q = qc* (antisymmetric with respect to the mirror 
plane) . 
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The pseudo spins r and a represent the signs of the amplitudes of the two relevant 
SALMs. For convenience, we shall henceforth call the latter ,r-SALM' and ,<t-SALM'. 
Previously || |J, we deduced the displacements corresponding to the r- and the a- 
SALMs for all BCCD atoms from experimental data [fT5|] . For a complete visualization 
of these SALM displacements it is sufficient to follow the ^-displacements of the nitrogen 
atoms. Consider the two nitrogen atoms 1 and 2 associated to a half cell. The major 
part of their ^-displacements is provided by the r-SALM contributions (equal for 1 and 
2). The (x-SALM contributions (of equal sizes but opposite signs for 1 and 2) account 
for the small relative y-displacement of 1 and 2. 

Taking the transformation properties of the SALMs into account, the form (1) of 
the pseudo spin Hamiltonian is obtained as well as expressions for the spontaneous 
polarization in terms of the pseudo spins || || . 

Experimental investigations of the pressure-temperature {p-T) phase diagram of 
BCCD [H], [18| revealed a dielectric anomaly in the modulated phase with wave 
vector k = |c*. It occurs along a line T$(p), which runs almost parallel to the boundary 
of the para phase and to the lines of constant spontaneous polarization Ps, x - At ambient 
pressure, Ts practically coincides with the upper boundary of the fourfold phase. The 
T^-anomaly exhibits the characteristics of a first order phase transition; it leaves the 
wave number unchanged. One proposal for the explanation of the T^-anomaly was to 



interpret it as a freezing of the soliton lattice [|T1|, i.e. to assume that the soliton lattice 
is free to move above Ts and is pinned below. 

Since the wave number k in BCCD is measured in units of the reciprocal lattice 
constant, whereas the wave number q in the DIS model is given in units of the reciprocal 
layer spacing, the k = |-phase in BCCD corresponds to the q = |-phase in the DIS 
formalism. For the latter the occurence of an internal transition was explicitely shown in 
sections |] and || above. It suggests itself to interpret the T^-anomaly as a manifestation 
of the internal transition observed in the DIS model. This hypothesis is supported by 
the following results: 

In the same way as the T^-line in BCCD, the X^Tine runs parallel to the 
boundary of the high-temperature phase and parallel to the lines of constant spontaneous 
polarization @ P S:X = T,t(-^Y[Pxt e t e+ i + P%stSt+i + P^ete + Pi (t e s e+1 - s e t e+1 )}. 

There is an interesting connection with experimental studies on the structure of the 
fourfold phase in BCCD: 



The displacement of atom /i in cell n can be written as (e.g. reference |19| ' 
< = J2 [ e m cos(2vrmk • (r n + ff)) + e^ 4 sin(27rmk ■ (r n + (f)) 



[ID 



where the zeroeth order harmonics are absorbed in the average positions p M of fi with 
respect to the origin of the cell. There are only odd harmonics due to the symmetry of 
the k = j-phase in BCCD. A comparison of the symmetries of the eigenvectors e^ 4 and 
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for odd m and of the r-SALM and the a-SALM shows that the sin and cos terms 
in (|TT|) correspond to the t- and s-profile respectively 0. 

Recent neutron scattering studies of the fourfold phase in BCCD at 100i^ 
and ambient pressure (see figure 0a), i.e. below T$ ~ 115K, revealed a strong third 
harmonic in the cos terms in (|TT]) . This is in perfect agreement with the structure of the 
low temperature modification in the DIS model as derived in sections |2] and [3[ X-ray 
techniques \T5\ showed a more sinusoidal modulation (T = 90K, ambient pressure), 
which can be well described as the high temperature modification in the DIS model, i.e. 
the one that is expected to be stable above T$. 

These seemingly contradictory experimental results can be explained if the influence 
of point defects on the phase diagram and recent experiments on intentionally X-ray- 
irradiated samples are taken into consideration. Exposure to X-rays - even to doses used 



normally in diffractometry [17[ - produces radiation damage, which is sufficient 



to cause an observable lattice expansion. The effect is comparable to doping with 
bromine ions. The resulting (positive) plastic strains shift the boundaries in the p-T- 
phase diagram to higher pressures/lower temperatures. It was experimentally observed, 



that this shift is stronger for the Tg-line than for the other phase boundaries [17, 18 



This leads to the situation described schematically in figure 0b: the irradiation defects 
induced by the X-ray measurements produce different shifts of the T^-line and the other 
boundaries and thus open up a window for the observation of the high temperature 
modification of the k = |-phase at ambient pressure. 

For a more quantitative comparison of the structures predicted by the symmetry- 
based DIS model with the atomic displacements observed either by X-ray or by neutron 
measurements, we derived the amplitudes Q\ and Ql in a decomposition of the 
structural modulation in the fourfold phase of BCCD in terms of the r- and a-SALMs 



(see figure §) from the experimental data given in references [19], [15|]. In the employed 
procedure ||, Q\ and QJ are proportional to ti and se in linear approximation. Q\ and 
Ql can be written in terms of first and third order harmonics as in equations (4) and (5). 
We obtained 03/01 = 0.15, 63/61 = —0.08 from the X-ray data |15| and 03/01 = 0.27, 
63/61 = 2.35 from the neutron results [19|. Similar values can be found in the DIS model 
for temperatures not too far above (e.g. 9 = 4.25: 03/01 = 0.18, 63/61 = —0.07) or below 
6*int (e.g. 9 = 2.76: 03/01 = 0.39, 63/61 = 2.13) respectively. It is noteworthy that even 
with the parameters (A = 0.05, j = j' = —1.0) which were rather arbitrarily chosen in 
this comparison and which probably do not represent the optimal choice (e.g. by the 
methods discussed in ref. ||) for BCCD, we find a good agreement between theoretical 
prediction and experimental observation. 
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5. Conclusions 

Recent theoretical developments in the field of modulated systems allow, in a physically 
well motivated way, to bridge the gap between pseudo spin model-parameters and 
-properties on the one hand, and experimental control parameters and physical 
characteristics of real systems on the other. 

An efficient pseudo spin model, which is well suited as a basis both for largely 
analytical calculations and for numerical simulations, is the DIS model. It can 
advantageously be used for the theoretical description of a large class of modulated 
materials. Since the model variables are explicitely related to local properties of the 
discrete crystal lattice, it allows explicit predictions e.g. of structures, symmetries, 
spontaneous polarizations, phase diagrams, and orders of transitions of modulated 
phases in real systems. 

A special feature of the DIS model is that it predicts phase transitions not 
only between modulated structures with different wave numbers, but also between 
modulations with the same wave number but different pseudo spin configurations 
('internal' phase transitions). The respective phases differ significantly in the amplitudes 
of the harmonics in a Fourier expansion of the pseudo spin profiles and in symmetry. 
Since remotely similar theoretical results derived for the ANNNI model were later shown 
to be artefacts of the employed mean field approximation, we verified the validity of our 
MFA results both by means of the mean field transfer matrix method and by Monte 
Carlo simulations. 

One indication of the suitability of the symmetry-based DIS model for the 
description even of detailed characteristics of complex modulated real systems is the 
proposed explanation of the dielectric anomaly T5 observed in the fourfold phase of 
BCCD: this anomaly and its characteristic properties can be well interpreted as a 
realization of the internal phase transitions of the DIS model. 

Furthermore, an application of the general procedure for translating model 
properties to characteristics of real, experimentally investigated systems allows a 
tentative explanation of discrepancies between experimental structure determinations 
of the fourfold phase of BCCD by neutrons on the one hand and by X-rays on 
the other. Taking the influence of X-ray induced defects (via the strain field) into 
account, we propose a mechanism which allows the high temperature modification of 
the fourfold phase (which exists in ideal crystals only at elevated pressure) to become 
stable already at ambient pressure and, thus, be visible in X-ray measurements. The 
two experimentally determined structures of the fourfold phase agree fairly well with 
our calculated data. 
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Figure captions 



Figure 1. Reduced temperature- interaction phase diagram of the DIS model with 
A = 0.1 and j = j' = —0.25 (see text). Shown are the paraelectric, the ferroelectric, 
and a few commensurate phases. Higher order commensurate and incommensurate 
phases are contained in the hatched areas. LP is the Lifshitz point, MP the multiphase 
point. 9 is the reduced temperature ^jj-, whereas k_ is given by = | (k — A). The 
phases are characterized by quotients q = m/n giving the modulation wave numbers 
in units of 27r/(pseudo spin spacing). 



Figure 2. Pseudo spin profiles for (a) the (12)-phase (q = g) and (b) the (2)-phase 
(q = |) as functions of the reduced temperature 9 with A = 0.05 and j = j' = — 1. 
Shown are the mean values t e (solid lines) and sg (dashed lines), for layers I from 1 to 
the modulation length, i.e. from 1 to 6 for case (a) and from 1 to 8 for case (b). 



Figure 3. Pseudo spin profiles tt (upper part) and sg (lower part) of the q = \- 
phase for reduced temperatures slightly below (9 = ksT/M = 2.967, left hand side) 
and slightly above (9 = 2.993, right hand side) the internal transition 9i n t- Only first 
(solid lines) and third (dashed lines) order harmonics in sine (cosine) are present in 
the Fourier expansion of tg (sg). 



Figure 4. Amplitudes of the first and third order harmonics in the pseudo spin profiles 
tg and si of the q = |-phase as functions of the reduced temperature 9 = ksT/M 
[see equations (4) and (5)]. The internal phase transition at 9i nt separates the low 
temperature modification with dominant third harmonic in s from the much less 
anharmonic high temperature modification. 



Figure 5. MFTM values for the reduced free energies F/M of the low- (dashed line) 
and the high-temperature (solid line) modifications of the q = | phase as function of 
the reduced temperature 9. Fixed parameters are A = 0.05 and j = f = —1. Solid 
line: low temperature modification; dashed line: high temperature modification. For 
better clarity, F/M - 9 + 4 is plotted. 



Figure 6. Specific heat of the q = i-phase obtained from Monte Carlo simulations of 
DIS systems with 30 x 30 x 6 spin sites. A = 0.1, k = -0.57, and j = j' = -0.5. The 
inset shows the calculated reduced energy as a function of 9. 



Figure 7. p-T phase diagram of BCCD near the k = 1/4-phase. Solid lines: 
boundaries of the fourfold phase; dot-dashed lines: Tg-anomaly. •: neutron scattering 
measurements [19]; *: X-ray diffractometry [15]. a) Undisturbed sample (data from 
reference [18]). At ambient pressure, T$ coincides with the upper phase boundary, 
b) Sample for which X-ray induced defects are (schematically) taken into account. 
The temperature decrease is larger for Tg than for the phase boundaries. X-ray 
measurements at * should establish the structure above T$ even at ambient pressure. 



Figure 8. Fourfold phase of BCCD: amplitudes Q[ and of the r- and cr-SALMs as 
determined by the corresponding decompositions of neutron [19] (left hand side) and X- 
ray [15] (right hand side) data. In a Fourier series expansion only first and third (solid 
and dashed lines respectively) order harmonics in sine (cosine) are present in the spatial 
modulation of Q\ (Q°). The neutron scattering results [19] show a much stronger third 
harmonic in the a-SALM-profile than the X-ray diffraction measurements [15]. 
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